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We present a preliminary calculation of Bk using domain-wall valence quarks and 2+1 flavors 
of improved staggered sea quarks. Both the size of the residual quark mass, which measures 
the amount of chiral symmetry breaking, and of the mixed meson splitting A m ; x , a measure of 
taste-symmetry breaking, show that discretization effects are under control in our mixed action 
lattice simulations. We show preliminary data for pseudoscalar meson masses, decay constants 
and Bk- We discuss general issues associated with the chiral extrapolation of lattice data, and, as 
an example, present a preliminary chiral and continuum extrapolation of f K . The quality of our 
data shows that the good chiral properties of domain-wall quarks, in combination with the light 
sea quark masses and multiple lattice spacings available with the MILC staggered configurations, 
will allow for a precise determination of Bk- 
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1. Lattice calculation of Bk with a mixed action 

The kaon B-parameter (Bk), which parameterizes the hadronic part of neutral kaon mixing, 
plays an important role in flavor physics phenomenology. When combined with an experimental 
measurement of Sk, B K constrains the apex of the Cabibbo-Kobayashi-Maskawa unitarity triangle. 
Because Ek is well-known, the dominant uncertainty in this constraint is that of B K . 1 It is likely that 
whatever new physics exists has additional CP-violating phases; these will manifest themselves as 
inconsistencies between measurements that are predicted to be identical within the Standard Model. 
Thus a precise determination of Bk will help constrain physics beyond the Standard Model. 

The determination of B K is an important goal of the lattice QCD community. Thus many calcu- 
lations have been done, each improving upon the previous one. The benchmark JLQCD calculation 
involved a thorough study of the quark mass and lattice spacing dependence, but did not include 
the effect of sea quark loops, resulting in an indeterminate quenching error [jl|]. The HPQCD col- 
laboration eliminated this uncertainty by using dynamical staggered fermions, but their result has 
a ~ 20% systematic error from neglected higher-order and wrong-taste operators in the lattice-to- 
continuum matching procedure The RBC and UKQCD Collaborations recently calculated Bk 
with dynamical domain-wall fermions and fully nonperturbative operator renormalization, but they 
have only a single lattice spacing and cannot yet perform a continuum extrapolation [^J. 

Our mixed action calculation combines domain-wall valence quarks and staggered sea quarks. 
We use the MILC 2+ 1 flavor improved staggered lattices which are publicly available with a large 
range of quark masses, lattice spacings, and volumes and allow for good control over the system- 
atic error from chiral and continuum extrapolation We use the domain- wall propagator code 
from the Chroma lattice QCD software package Because domain- wall quarks do not carry 
taste quantum numbers, the Bk lattice operator mixes only with other operators of incorrect chi- 
rality, making the chiral extrapolation more continuum-like than in the purely staggered case, and 
allowing for nonperturbative operator matching using the Rome-Southampton method [|6|]. Thus 
the mixed action method combines the advantages of staggered and domain-wall fermions without 
suffering from their primary disadvantages and is well-suited to the calculation of Bk- 

2. Discretization effects in mixed action lattice simulations 

Each flavor of staggered quark comes in four identical species, or "tastes", that are related by 
an SI/ (4) symmetry in the continuum [{7|]. The staggered taste symmetry, however, is broken at 
nonzero lattice spacing, causing the masses of the sixteen tastes of sea-sea mesons in mixed action 
simulations to be split according to their SO(4) taste representation at @{a-) ||]: 



These taste-splittings generically enter staggered #PT expressions for masses and weak matrix 
elements. They have already been measured, however, for the MILC lattices in Ref. [Q], and we 
use these values in the chiral and continuum extrapolation of our lattice data. 

'While Bk is currently the dominant source of uncertainty, the importance of other quantities such as the perturbative 
Inami-Lim functions and [V c j| will increase as the precision on Bk improves. 




}i(m s + m s i) +a 2 A t . 



(2.1) 



2 



Bfrfrom unquenched mixed action lattice QCD 



Ruth S. Van de Water 
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Figure 1: Chiral extrapolation of m res on the coarse MILC lattices. Because there is no notion of full QCD at 
nonzero lattice spacing in the mixed action theory, the curve shows the extrapolation/interpolation for points 
where the domain-wall pion mass is tuned to equal the lightest (taste pseudoscalar) staggered pion mass. For 
comparison, we show the determination of m re s by the LHP collaboration, which uses this tuning Jic[], 

Domain-wall quarks receive an additive contribution to their mass from explicit chiral symme- 
try breaking, the size of which is controlled by the length of the fifth dimension [^]. Consequently, 
to lowest order in %PT, the masses of valence- valence mesons in the mixed action theory are 

m vv , = /I (ray +my +2m res ), (2.2) 

where m res is the residual quark mass. In our simulations the length of the fifth dimension is 
L$ = 16. With this choice, we determine a preliminary value of r\m res = 0.00446(9) (45), around 
3 MeV, in the chiral limit; the extrapolation is shown in Figure [l| Because m res is a quarter the size 
of our lightest valence quark mass, and half that of RBC/UKQCD (in physical units), this indicates 
that chiral symmetry breaking is under control in our mixed action lattice simulations. 

Because the mixed action lattice theory has new four-fermion operators, the chiral effective 
theory can have new low-energy constants. It turns out, however, that the mixed action chiral 



Lagrangian has only one new constant at lowest order fjl 1|] . This coefficient, due to taste-symmetry 
breaking from the sea quarks, produces an ff(a 2 ) shift to the mixed valence-sea meson mass- 
squared: 

ra V5 = /I (ray + m res + m s ) + a 2 A mix . (2.3) 

In order to determine A m ; x for our lattice simulations, we construct a linear combination of valence- 
valence and valence-sea squared meson masses, 

^ {2m\ y - m\ s ) = a 2 A mix + 2/lra s , (2.4) 

and perform a linear fit versus the staggered quark mass; this fit is shown in Figure ^[ We obtain 
r^a 2 A m i x = 0.206(16)(21), where the systematic error is preliminary. This is close to the staggered 
sea splitting Aa. It is also consistent with the independent determination of A^ by Orginos & 



Walker-Loud [12]. Because we have data at two lattice spacings, we can also determine the scaling 
behavior of A m i x . We find that it scales correctly for the MILC lattices, as Gia^a 1 ). 
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circles: a=0.12fm; squares: a=0.09fm 




Figure 2: Determination of A m j x on the MILC coarse (circles) and fine (squares) lattices. The filled symbols 
are our data points. For comparison, the open symbols show the four MILC staggered taste splittings. 



3. Preliminary data and analysis 

We have generated quark propagators on the MILC coarse (a « 0.12 fm) and fine (a « 0.09 
fm) ensembles with valence quark masses from ~ m s /S - m s . In order to keep finite volume effects 
under control, we restrict the combination am n L > 4; our lightest pion is ~ 280 MeV. Unlike in 
other mixed staggered sea, domain-wall valence simulations, we do not tune our valence-valence 
pion masses to any particular value. Because there is no notion of full QCD at nonzero lattice 
spacing in the mixed action theory, we generate many partially quenched data points, and we use 
the appropriate partially quenched, mixed action #PT to extrapolate our lattice data. Our analysis 
is still in progress, and in this section we present only preliminary data and chiral fits. 

A necessary prerequisite for the calculation of Bk is the determination of the pseudoscalar 
decay constants f % and ■ Because they are known quantities, they provide a test of the lattice 
methodology, and allow one to demonstrate an understanding of and control over systematic errors. 
In order to suppress contamination from pions circling the lattice in the time direction, we use 
symmetric and antisymmetric linear combinations of periodic and antiperiodic quark propagators 
in our correlation functions. We then extract the decay constant using the axial Ward identity: 

_ A WP V2(m x + m y + 2m res ) ^ ^ 

where Aww, A^p, and m n come from a simultaneous fit of wall-wall and wall-point correlators. 

Figure || shows preliminary chiral and continuum extrapolations of the "pion" mass and decay 
constant data (i.e., m x = m y ) using the approach of the MILC collaboration described in Ref. tjl3|]. 
They are each simultaneous fits to the coarse and fine data using the next-to-leading order (NLO) 
partially quenched, mixed action #PT expression for m n or f„, plus additional NNLO analytic 
terms. In order to achieve good correlated fits without needing NNLO logarithm terms (the full set 
of which are not known for the mixed action theory), we omit valence points with m n > 500 MeV 
and all points on the heaviest {mi/m s = 0.4) coarse ensemble. We also correct the data by the 
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Figure 3: "Pion" mass-squared and decay constant versus average quark mass. The circles are coarse data 
points and the squares are fine data points. The cyan band is the continuum full QCD curve (with statistical 
errors) that results from fitting the lattice data to the appropriate mixed action ^PT expression; fit (a) has a 
CL = 0.34 and fit (b) has a CL = 0.83. Note that the band in plot (b) contains the experimental value of f n . 



known 1-loop finite volume effects. Note that, with these judicious cuts, we have 30 data points for 
only 7 fit parameters, and do not need to introduce Bayesian constraints. 

Although the fitting procedure described above successfully reproduces the experimental value 
of f n , it is not the only reasonable possibility, and various physically motivated approaches are 
used in the literature. One can, in principle, continue to remove heavy masses until unadulterated 
NLO #PT gives an acceptable fit. Typical lattice simulations, however, only have one or two 
data points in this regime, making this method largely impractical. Furthermore, one cannot do 
this for quantities involving m strange . One can instead include heavier data points at the cost of 
introducing terms of higher-order in #PT. Unfortunately, however, the NNLO chiral logarithms are 
known only for the pseudoscalar masses and decay constants and only in the continuum [JlJ]. It 
is likely that the inclusion of only NNLO analytic terms does not significantly impact extrapolated 
results for physical quantities [15], but, of course, it would be preferable to include the full NNLO 
expression. Finally, use of SU (2) #PT provides a valuable crosscheck on the determination of m n 
and f K , but cannot be used to determine quantities involving valence strange quarks (such as f% 
and Bk ) without further assumptions. One can either match an SU (2) #PT fit of light quark masses 
onto another fit of data near m strange , which requires an additional fit ansatz, or one can treat the 
strange quark as a heavy quark within SU (2) heavy meson #PT, but this relies on the assumption 
that m str ange is heavy compared to Aqcd- Ultimately, the degree to which these fit choices change 
results for extrapolated physical quantities must be reflected in the systematic errors. 

Although we do not yet have sub-percent statistical error bars like the MILC Collaboration, 
we can still begin to distinguish between fit ansatze that work and those that do not, and hence draw 
interesting conclusions about the application of %YY to lattice QCD. In order to assess whether or 
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Figure 4: Lattice Bk versus light quark mass in "kaon". 



not we have sufficient data in the chiral regime to allow the use of #PT, we must first look for chiral 
logarithms in the pseudoscalar data. We observe the appropriate curvature in the quantity m\jm q , 
which is expected to have large partially quenched chiral logarithms. We also find that, in order to 
get a fit of the pion decay constant data with an acceptable confidence level, we must include the 
NLO chiral logarithms. Thus we have both qualitative and numerical evidence for the presence of 
chiral logarithms in our data. Not only must we include chiral logarithms for good fits, but in fact 
we need the correct mixed action chiral logarithms. We can no longer fit the decay constant data 
if we set the taste-breaking mass-splittings A/, A m i x — > 0. Thus we conclude that it is important 
to use the %PT expression appropriate to the specific lattice action used in simulations. 2 We also 
implicitly observe finite-volume effects in our data. Despite our reasonably conservative choice 
to keep the quantity m n L > 4, the 1-loop finite-volume corrections to our m\ data are as large as 
~ 3% due to enhancements from partial quenching, and we cannot get an acceptable fit to our pion 
mass data without their inclusion. Thus we are currently performing a finite-volume study on one 
of the coarse ensembles. Finally, our decay constant fit clearly resolves a large &{a 2 ) analytic term 
that causes the coarse data to lie ~ 10% above the continuum full QCD curve. Because generic 
discretization effects will likely be of this size for all simulations not using &(asa 2 ) improved 
lattice actions, this underscores the importance of having multiple lattice spacings. 

We calculate the K° — K° matrix element using Coulomb gauge-fixed wall-source propagators, 
and then determine Bk with the following ratio: 



B K 



V 



*-wpw 



if source j t j tsink ) 



8/3 Cm p{t source -i^Cyj pit s i n k 



(3.2) 



where, in the upper labels, "P" indicates the pseudoscalar operator, u ff" indicates the 4-quark 
operator, and "A" indicates the axial operator. Figure ^| shows the resulting preliminary data. 



2 Another quantity where inclusion of the correct lattice discretization effects is necessary is the isovector scalar 
correlator. We cannot describe the behavior of our ao correlator data jl6| if we set A m ; x — > in the mixed action ^PT 



expression for the disconnected bubble contribution Jl7[ ] . 
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4. Summary and future plans 

Br must be known to an accuracy of 5% or better in order to have phenomenological impact. 
We have already calculated the expression for Bk in mixed action ^PT and shown that, at NLO, it 



has only one more low-energy constant than in the continuum []18[]. We have determined the sizes 
of m res and A m ; x on the MILC lattices and shown that both the amount of chiral symmetry breaking 
and taste-symmetry breaking are under control. We have generated correlation functions necessary 
for the calculation of pseudoscalar masses, decay constants and Bk at several sea quark masses and 
lattice spacings, and have begun chiral and continuum extrapolations. In the near future we plan 
to finish the nonperturbative renormalization and determine the domain-wall valence strange quark 
mass. Our preliminary results are promising. The errors in our Bk data range from 0.5-2% and 
will continue to improve, and our preliminary determination of f„ has a 2.5% statistical uncertainty. 
This suggests that the mixed action method will allow a precise determination of Bk- 
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